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Introduction. 

In [4] Bailey et al (among other things) study the Bessel moments 



Cm.k = y ^ Ko{x)"^dx 



Here Ko{x) is a certain Bessel function that conveniently can be defined by 

oo 







This leads to another representation (in Ising theory) 



oo oo 



Cm,k — „„,, 



k\ f f dxi...dxr, 



2™ 7 "J (cosh(a;i) + ... + cosh(a;™))'=+i 



(historically it was the other way around). 

In J. Borwein-Salvy [5] recursion formulas for the Cm,k are derived (m fixed). 
In the first section these recursions are studied in more detail. E.g. if we define 

_ 16" 

«ri — — ^7*^4,271+1 

we find an Apery-like recursion (compare [3]) and recognize formulas from [1] 
and [3]. Similar transformations of C5^2n+i lead to a 4-th order differential 
equations whose mirror at a; = oo is a Calabi-Yau equation found by Verrill 
(#34 in the "big table" [2]). This is also the case with C6.2ri+i where the 
differential equation at oo is of order 5 (also found by Verrill) with a Calabi-Yau 
puUback of order 4 (#130 in [2]). 

There is an infinite sequence of differential equations of Verrill where the 
coefficients are 



l/l ijji . 



In [6] she gives a rather complicated formula for computing the recursion of 
In the second part we simplify this essentially using ideas in J. Borwein- 
Salvy [5]. 

In the last section we prove the 

Main Theorem For m ^ 3 we have 



oo 

Y — 



4"n!2 



and 

■CO 

n=0 

satisfy the same Picard-Fuchs differential equation of order m+ = to/2 if m 
is even and = (m+ 1)/2 if to is odd. This equation is easily found by a Maple 
program. 

There is a simplified version of this result for Bessel fans: 
The differential equation satisfied by 

oo 

also has the solution 

w = x-'^Io{x-'^r 

This depends on the identity 

I. Some examples. 
Four Bessel Functions 

On p. 13 in [4] Bailey et al define 

oo 


where Kq is a Bessel function. In [5] the following recursion is derived 

64(fc + 3)c4,fe+4 - 4(A; + 2){hk'^ + 20k + 23)c4,fe+2 + {k + ifc^^k = 
We make the substitution 

16" 

dn = — 79-C4,2n+l 

nr 

and get the recursion 

(n + 2fdn+2 - 2{2n + 3){5n^ + 15n + 12)(i„+i + 64(n + ifdn = 

Then 

OO 

satisfies the differential equation where = x-^ 

- 2x{20 + 1){59^ + 50 + 2) + 64x2(6* + 1)^ 



2 



which we recognize as equation (a) in [1]. Then 

A. 



^ fnV f2k\ f2n-2k 



f-^ \kj \ k J \ n — k 

fe=0 \ / \ / \ 

satisfies the recursion with initial values A_i = 0, Aq = 1. Let -B„ be the 
solution with Bq = 0, i?i = 1. Then we have 
Theorem. We have 

dn = g^nC(3) - 3B„ 

Proof. In [4] we find 04^1 = |C(3) and 04^3 = ^C(3) — jq giving 
do = |C(3) and di = |C(3) — 3 Then we use the recursion. 

We want to find the asymptotic behaviour of An and d„ as n ^ 00. Making 
the Ansatz 

An = Cn^y 

in the recursion we find A = 16 or A = 4 and 6 = — |- Numerical experiments 
suggest 

16" 

An ~ 0.36- 



and 

This gives 
which proves 



dn ~ 0.7 



n3/2 
4" 



3/2 



24^^ ^ A„ 4" 



Remark. The differential equation 

- 2x(20 + 1)(56»^ + 50 + 2)+ 6Ax^{e + if 
is self dual at infinity and the coefficients can be written (H.Verrill, [6]) 



z+]+k+l=n •' 



Five Bessel functions. 

Consider 



C5,2„+i = j x'^'^+^Koixfdx 
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Then using the ideas of [5] we find the recursion 

225c5,„+6 - (259n^ + 1554n + 2435)c5,„+4 

+(35n^ + 280n^ + 882n^ + 1288n + 731)c5,n+2 - (n + lfc5,n = 
Make the substitution 

dn — pr'^5,2n+l 

which gives the recursion 

n^{n - l)2fi„ = 4(n - f )2(259n2 _ 5l8n + 285)d„_i 
-3600(35n'^ - 210n^ + ASSn^ - 504n + 201)rf„_2 + 3240000(n - 2)'*rf„_3 

Let An be the solution of the recursion with initial values = 1, = 0, 
A2 = 0. Similarly let iJ„ and C„ be solutions with Bq = 0, Bi = 1, B2 = 0, 
Co = 0, Ci = 0, C2 = 1 respectively. Then 

dn = Ans + 225Bnt + C„(6750 - 4500s + 64125t) 

where ,s = C5_i and t = 05,3 . We also use the conjectured value of 05,5 = 
— ||.s + Unfortunately we still do not know the exact values of s and 

t. Maybe they arc related to the Apery limits of ^ and ^ 
A related Calabi-Yau equation. 
With 9 = x-^ the differential equation satisfied by 

00 

y = Yl 

is 

6»^(6»-l)2-4x6'^(2596»^+26)+3600a;^(356»'*+706i^+636'^+286i+5)-3240000a;^(6i+l)^ 

The last factor cointains (6 + 1)'^ which suggests that transforming the equation 
to a: = 00 could give a Calabi-Yau equation. This is indeed the case: The 
substitutions 9 — > —9—1 and x — > 900a;^^ give 

9'^ - x{359^ + 709^ + 636*^ + 286» + 5) 

W{9 + if {2599^ + 5186 + 285) - 225a;^(6' + if {9 + 2f, 

an equation found by Helena Verrill [6] . It has #34 in the big table [2] and has 
the analytic solution 

00 

y = ^ a„a;" 

n=0 
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where 

Six Bessel functions. 

Consider 



C6,fe = y x^KQ{xfdx 



As above we have 

2304(fc + 4)c6,fe+6 - 16(A; + 3)(49A;^ + 294A; + 500)c6,fc+4 

+8(fc + 2)(7fc^ + 56A;^ + 182fc^ + 280fc + 171)c6,fe+2 - (fc + l)^C6,fe = 
With the substitution 

48^" 

rfn = 



we have the recursion 

(2n + 5)(n + 3)^(n + 2)2d„+3 - 32(n + 2)^(196n^ + 784n + 843)d„+2 
+64 • 48^(271 + 3) + 84n3 + IQen^ + 210n + 87)d„+i - 128 • 48^(71 + l)^d„ = 

Consider the three solutions An,Bn, C„ with initial values 

Ao = 1, Ai = 0, = 

Bo = 0, = 1, = 
Co = 0,Ci =0,^2 = 1 

respectively. Let C6,i = s, C6,3 = t . Then C6,5 = ^ — ggs + f|t is conjectured. 
Then we have 

dn = Ans + 2304B„i + C„ (138240 - 36864s + 1566720i) 

A related Calabi-Yau equation. 

Let 

OO 

y = 5Z 

n=0 

Then y satisfies the differential equation 

e^{0 - if {29 - 1) - 32x6'^ (1966*^ + 59) 
+64 • A8'^x'^{2e + 1)(146''^ + 286*^ + 286*2 + 146* + 3) - 128 • 48V(6» + 1)^ 



We find the mirror equation at a; = oo via the substitution 6 — > —6 — 1 and 

X — > 96^a;~^ 

- 2x{20 + 1){U0'^ + 286^ + 280^ + 146» + 3) 

+Ax^{e + 1)^(1966'^ + 39261 + 255) - 1152x^(6 + if {6 + if {26 + 3) 

This wc recognize as #130 in the big table. It was found by H.Verrill [6]. The 
coefHcients are 

A = V ( 1^ 
, 4^ HljlkUlmlsl' 

i+j+k+l+m+s=n 

Seven Bessel functions. 

Let 

, 105^" 

Then 

OO 

satisfies 

6^{6 - if {6 - 2f - 8x6'^{6 - 1)^(64586»^ - 64586* + 2589) 

+48 • 105^a;^6i^(6586»^ + 39661^ + 17) 
-64 • 105^x^(846'*^ + 2526'^ + 3786'^ + 3366^^ + 18061^ + 546' + 7) 
+256- 105V(6i + 1)'^ 
The transformation to infinity by 6 — > —6 — 1 and x — > 210^a;~^ gives 

6^ - x{8A6^ + 2526'^ + 37861* + 33661^ + 1806'^ + 5461 + 7) 

3x^(6' + 1)2(6586'* + 26326'-' + AMAO^ + 34246' + 1071) 
-2x^(6' + 1)2(61 + 2)2(64586'2 + 193746' + 15505) 
+1O5V(0 + 1)2(0 + 2)2(0 + 3)2 

with solution 

oo 

where 



, , i\j\k\l\m\p\s\' 

i+3+k+l+m+p+s=n 



II. Sums of squares of generalized binomial coefficients. 
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In [6] Verrill has given a rather comphcated formula for the recursion of 
We will instead consider 

X ^ 1 



»i+»2 + -..»fc=n 

Consider 

OO T 

Ex-' 

Then y satisfies the differential equation 

-X 

Actually 

y{x) = /o(4v^) 

Then 

OO 

n=0 

Using Lemma 3 in J.Borwein and Salvy [5] we find the following Maple program 
for computing the difi^erential equation for w for all m. 

S:=proc(m) local M,k; M(0):=1; M(l):=t; for k to m do 

M [k+ 1] :=x*difr(M [k] ,x) +M [k] *t-k* (m-k+ 1 ) *x*M [k-2] ; od; 
series(expand(M[m+l] ,x=0,infinity) ; end; 

Let m+ = m/2 if m is even and m+ = (m + l)/2 if m is odd. Then 
write 

Then the diflEerential equation satisfied by 



n=0 n=0 

is given by 



j=0 s=0 
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III. Proof of the Main Theorem. 

The Bessel function Ko{x) satisfies the differential equation Tm{x,6) given 
by the Maple program 

e^-x^ 

Using Lemma 3 in Borwcin-Salvy [5] again wc obtain the differential equation 
Tm{x,9) satisfied by Ko{x)"^ given by the Maple program 

T:=proc(m) local L,k; L(0):=1; L(l):=t; for k to m do 
L[k+l]:=x*difr(L[k],x)+L[k]*t-k*(m-k+l)*x*L[k-2];od; 
series(expand(L[m+l] ,x=0,infinity) ; end; 



The crucial part of the proof is the following 
Lemma. We have 

Mfc(x,e) = 2-('=+i)Lfc(2v^, 20) 
Proof: We use induction on k. Assume 

Mk-i = 2-'=Lfe_l(2^/i,2^) and Mk = 2-('=+l)Lfc(2^/i, 2^) 

Then 

dMk 

Mfc+i = x^^ + MkO - xk{m -k + l)Mk-i 
= x2-('=+i) ^Lfc(2v^, 26l)+2-('=+i)Lfe(2v^, 26i)6'-x2-'=A;(m-A;+l)Lfe-i(2v/5, 20) 
= 2-(fe+i)a;^— ^Lfc(2^/i, 20)+2-^''+^^Lk{2^, 20)20-{2^)H-^>'+^h{m-k+l)Lk-i{2^, 20) 



= 2-('=+2)Xfc+i(2v^,2^) 



The rest of the proof is merely book-keeping. Recall that 

T^{x,0)=J2^^'Pj{O) 

3=0 

annihilates Kq{x)"^. Then by the Maple program following Example 5 in [5] 

we find the recursion for c„i,k by substituting 6 — > — fc — 1 — 2j in Pj{0). 
Since k = 2n+l we get — > — 2(n + 1 + j) Then with 



8 



we get the following recursion for 



m-|- 

^ n'^in + lf...{n + jf4"'+-^Pj(-2{n + 1 - j))N^ = 

where N f(n) — f{n). Converting to the differential equation for y — rf„a;" 
we make the substitution n — > 9 — j and — > a;™+ in the coefficient of 

m-|- 

Y^x'^+'^e^ie - 1)2. ..(61 - jfA"'+-^Pj{-2{9 + 1)) 

To get the differential equation at oo we make the substitution 6 — > —6 — 1 
and X — *■ and we get 

m-|_ m+ 

which is the differential equation satisfied by 

oo 
n=0 
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